P=[1 2;2 3];
D=[1 0; 0 3];
A4=P*D.^4*inv(P)